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Abstract 

We construct an extended Hubbard model with open boundaries from a i2-matrix 
based on the U q [Osp(2\2)] superalgebra. We study the reflection equation and find two 
classes of diagonal solutions. The corresponding one-dimensional open Hamiltonians are 
diagonalized by means of the Bethe ansatz approach. 
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1 Introduction 



In recent years there has been substantial research devoted to the study of (l+l)-dimensional 
solvable lattice models with integrable boundary conditions. In one-dimensional theory with 
factorized scattering, the boundary effects can be described by scattering matrices satisfying 
the so-called reflection equation A systematic approach to construct such models has been 
developed by Sklyanin || who has investigated the six vertex model with boundary fields. 
Subsequently, this scheme has been generalized to handle a rather general class of models 
based on Lie algebras [[3], |J (see also ref. ||). By now, variants of this method have been 
extensively used in the literature to study various integrable quantum field theories and lattice 
statistical mechanics models with open boundaries. See refs.|| |jj just to mention few examples. 

Of particular interest are supersymmetric generalizations of the Hubbard model H |9|, |I~0|j 
due to their possible relevance in describing strongly correlated electron systems. These models 
are often derived from supersymmetric solutions of the Yang-Baxter equation invariant by the 



gl(2\l) and Osp(2\2) Lie superalgebras [II], [Hj]. A successful example is the supersymmetric 
free-parameter model with open boundaries constructed from the .R-matrix associated to the 
four dimensional representation of gl(2\l) ||13|| . The purpose of this paper is to perform similar 
task for another Osp(2\2) .R-matrix solution found some time ago by Deguchi et al fTI| . This 
latter model, however, appears to be specially interesting, since its Hamiltonian provide a 



lattice regularization of a relevant integrable double sine-Gordon model []14| . This continuum 
field theory with open boundary is known to describe tunneling effects in quantum wires [T5 



Therefore, the Bethe ansatz results for open boundaries we shall derive here could be a useful 
non-perturbative tool to investigate this condensed matter system as well. 

This paper is organized as follows. The next section is concerned with diagonal solutions 
of the reflection equations Jl|, § associated with a particular Osp(2\2) i?-matrix. We found 
two families of one parameter solutions leading to four possible choices of boundary conditions 
for integrable open chain Hamiltonians. In section 3 we formulate their Bethe ansatz solutions 
in terms of the coordinate Bethe ansatz approach. Section 4 is reserved for our conclusions 
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and final remarks. For completeness, in appendix A we present another possible i?-matrix 
embedding and discuss its boundary behaviour. 



2 The Osp(2\2) R- matrix and related K- matrices 

The "bulk" part of (l + l)-dimensional integrable system with a boundary is governed by the 
two-particle scattering matrix R(X) satisfying the Yang-Baxter equation. The spectral param- 
eter A plays the role of the difference of the particles pseudomomenta. The boundary effects 
are described in terms of two boundary K±(\) matrices that characterize the scattering in- 
teractions at the boundary |], 0. The compatibility of the reflections and particle scatterings 
leads us to an algebraic condition, the so-called reflection equation 0, Q 

Rn(X ~n)K (A)#2i(A + fi)K (//) =K 0)-Ri 2 (A + fi) K (A)P 2 i(A - //) (1) 

Here we will also require that i?-matrix i2i 2 (A) satisfies, besides the standard properties 
of regularity and unitarity, certain extra relations denominated PT and crossing symmetries, 
namely 

PT - symmetry : Pi 2 Pi 2 (A)P 12 = R\\ t2 (A) (2) 

crossing — symmetry : Pi 2 (A) = ^[ - — - V R%(— A — (3) 

p(-\-ri) 

where Pi 2 is the exchange operator, tk denotes the transpose in the space k, rj is the cross- 
ing parameter, V is related to the crossing matrix by M = V l V and p(A) is a convenient 
normalization function. 

When these properties are fulfilled one can follow the scheme devised by Mezincescu and 
Nepomechie [Q. One of the boundary matrices, say the i^-(A) matrix, is obtained by solv- 
ing the reflection equation (1). The other K + (X) matrix is automatically determined by the 
following isomorphism 

K+(\) = K t _(-\-7i)M (4) 

1 We are assuming the ordinary non-graded boundary framework, since we shall look only for diagonal 
solutions. For an application of the graded formalism to non-diagonal solutions, see for instance ref. |19[. 



Before proceeding, we remark that the original Osp(2\2) i?-matrix given in ref. JTl is not 
PT- symmetric. However, this difficulty is easily circumvent by applying a convenient unitary 
transformation. It turns out that the PT-symmetric Osp(2\2) i?-matrix can be presented in 
the following form 



R(X) 



q e*-q, x 



a=l 



q a l - eV 



1 — e A g 2 



a, (3 = 1 
a^P,p' 



1-q 2 



1 — e x q 2 



a, /3 = 1 a, 13 = 1 

a<(3,aj£(3' a>(3,a^/3' 



l X - 1 



1 - eV a =l 



a af3 (\)E af3 <g> E t 



' a 1 

a [3 



(5) 



where q is the deformation parameter, E a p are the elementary 4x4 matrices and we set 
a = 5 — a. The functions a a p(\) are 



<? Q e A + 1 



-i _ii-<? 2 



e A - 1 a,B 



Q /3 l + e A 

and the parameters q a and e Q are defined by 



1°C0' 



a = (3, 
a < (3 
a> (3 



(6) 





(7) 



It is not difficult to verify that this i?-matrix satisfies the properties (2,3), where 

/ A 



r) = in, p(A) 



1 



1 — e x q 2 ' 



V 



10 
10 
V -i / 



Now we have the basic ingredients to built up integrable boundary models associated with 
the i?-matrix (5). In this paper we are interested to look for diagonal solutions of the reflection 
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(9) 



equation, which will play the role of boundary fields in the context of spin chains. Substituting 
the ansatz 

/M«(A,£_) \ 

£«(A,f_) 

C«(A,f_) 

\ D«\\,{;_)J 

in equation (1), and after some algebra we find two classes of diagonal solutions given by 

A«(A,£_) = (e A + e«-)(e A -e«-g- 2 ) 

flW(A,e_) = (e- A + e«-)( e A - e «-g- 2 ) 

C«(A,£_) = ( e - A + e «-)(e A -e«-g- 2 ) 

£>«(A,^_) = ( e - A + e «-)(e- A - e «-g- 2 ) 



(10) 



and 



A( 2 )(A,e-) = (e A + e«-)(e A -e«-g- 2 )(e A -e^-g 2 ) 

5 (2) (A,£_) = ( e - A + e «-)(e A -e«-g- 2 )(e A -e-«-g 2 ) 

C( 2 )(A,e-) = (e A + e^)(e A -e^g- 2 )(e- A - e ^-g 2 ) 

D {2 \\,£_) = (e x + e^-)(e' x -e^-q~ 2 )(e~ x -e^-q 2 



;n) 



The corresponding fT^(A,£ + ) matrices are easily derived from the isomorphism (4) with 
the help of crossing property (8). However, it is convenient to introduce suitable K+)(\,£ + ) 
matrices to make the interactions at the boundaries be as much as possible symmetric. With 
this in mind, we choose the following expression for these matrices 



iif (A,^) = K w (-\ - in, in + 21og[g] - £+)'M 



'(0 



(12) 



where £± are two independent parameters that characterize the interactions at the right and 
left ends of the open chain, respectively. 

Equipped with the refection matrices K±(\,£±), an integrable model with open boundary 
condition can be obtained through the double-row transfer matrix t(X) formulated by Sklyanin 



E 



n (m) o (I) 

k\ (A)r(A) k! {\)t-H 



-A) 



m 



1,2 



(13) 



where T(A) = i? a ^(A) ■ • -i2 a i is the standard monodromy matrix of the corresponding closed 
chain with L sites. We note that the fact K+. (A, £+) can be taken from either if_ (A,£_) or 
ICi (A,£_) gives us four possible choices of boundary transfer matrices. In next section we 
will present the Bethe ansatz solution of the corresponding four one-dimensional Hamiltonians 
with open boundaries. 



3 Integrable open boundary Hamiltonians 



In order to obtain solvable Hamiltonians with open boundaries we have to expand the double- 
row transfer matrix t^' m \X) up to the second order in the spectral parameter A fI6 |. This is 
because the relevant term in the first order expansion, which is proportional to the trace of 
.K+^O), turns out to be zero for both classes of i^-matrices of section 2. Following ref. it 
can be shown that the corresponding Hamiltonian commuting with t"' m \X) is given by 



L-l 

i T r 

2 a 



' 2dX ' £?M 



Tr r , 



d 



a A 



-I Tr 



(0)F 



(14) 



where 



^■ +1 = p Jij+1 ^E Jii+1 (o), i? 12 (o) = CP12, ^ (m) 



Tr„ -A-i""(0) + ^Tr a K^(0)H La 

(15) 

We recall that in the derivation of expression (14) it is implicitly required that K^(X, £_) 
has been normalized by a scalar function such that if_ (0, £_) = Id. For both cases (10,11), 
the only non-trivial contributions to the boundaries come from the first two terms of equation 
(14). The remaining ones are proportional to the identity, playing the role of irrelevant additive 
constants. Considering the physical motivations of the introduction, one would like to express 
these open Hamiltonians in terms of fermionic creation and annihilation fields cj CT and Cj jCr 
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acting on the site j and carrying spin index a — ±. After some algebra, the final expression 
for the Hamiltonians can be written as follows 

H ii,m) = J2 Hj j+1 + U(C-)n 1+ni _ + U{i + )n L+ n L _ + £ ^n la + ^n Lc (16) 

3=1 <r=± 

where n^„ = Cj a Cj t(T is the number operator for electrons with spin a on site j. The expression 
for the bulk part is 

H 3,3+i = Yl [ c L c i+L- + h - c ] t 1 ~ nj-'O- + aV i) ~ n 3+h-^ 1 ~ aV i)\ 



a=± 



+ iVi \ n io - n i+lo-} 

<T=± 



+ V 2 [c) >+ c)_C j+1 ,-Cj+i )+ - c) >+ c) +l _c j+1 , + Cj _ + /i.e. 

where the couplings Vi and V 2 are determined in terms of the parameter q = exp(ry) by 

VI = sin(7), V 2 = cos(7). (1 



+ n j+ i i+ n j+ i _ + n j)+ n j+ i _ + n j5 _n i+ i i+ - ^ (n ja + n j+1<T ) 

a=± 



(17) 



Turning now to the boundary interactions we found that the on-site Coulomb coupling 
U{£±) is given by 

U(t±) = -t- 



sin(27) 



sinh(e±/2 - i 7 ) cosh(e±/2) 
while the boundary chemical potentials are 



and 



(m) 



,.(1) _ ..(1) _ ivl 

P1+ — An- — lv l 
/4+ = Wx 



e -«_/2+i 7 
sinh(g_/2-i7) 
C-/2+i7 (2) 



sinh(^_/2-i7) 



, til = -iVi 



cosh(£_/2) 



(i) (i) . T/ e e+/2- t 7 

,( 2 ) 



sinh(£ + /2-j 7 ) 
£ + /2-«7 _ (2) 



..W _ ? -T/— £— — — // l ; — 7'V,- 



-S+/2 



(19) 



(20) 



(21) 



L cosh(£ + /2) 

Our next task is to diagonalize the Hamiltonian (16) by the coordinate Bethe ansatz for- 
malism. The number of electrons N° with spin a are conserved quantities and they label 
the possible disjoint sectors of the Hilbert space. For a sector of a given number of particles 
N e = + N~ the Bethe wave function assumes the following form 



X Qj > °J P 



3=1 



,-,^ e 4 Ql -i., |0> (22) 
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where |0) denotes a reference state containing none particles, 1 < Xq 1 < xq 2 < • • • < xq n < L 
indicate the positions of the electrons, P is the sum over all the permutations of the momenta 
(Pi • ■ ■ Pjv e ) and the symbol sgn denotes the sign of the permutation. For configurations such 
that \xq i — xq^ > 2 the Hamiltonian (16) behaves as a free-theory and the solution of the 
eigenvalue problem \ty = E^ l ' m \L) \$) is 

E^ m \L) = 522cos(%) (23) 
i=i 

up to some additive constants proportional to the number of electrons N^. It is standard in 
Bethe ansatz approach that configurations in which the electrons are nearest neighbors, at the 
same site or even at the boundaries impose constraints on the amplitudes of the wave function. 



For previous similar computations to other models with boundary, see refs.[17, 18]. We found 



that such consistency condition on the "bulk" provide us to the following relation 

^■■■o-j, CTi - (' ■ ■ ) kj, ki, ■ ■ •) = Sij(ki, fcj)A...,j i)0 y..(- • • , ki, kj, ■ ■ •) (24) 

while the reflection at the left and right ends of the chain gives us 

-/.-,•••) = S^P^A^...^...) (25) 
A.., ai (•••,-%) = S I {k ji P^)A.^X'" i k j ) (26) 

The two-body 5*-matrix Sij(ki,kj) connects the scattering amplitudes between the states 
{(fcj,<7j); (kj,<7j)} and {(kj,a'j); (fcj,cr-)} and its non-null elements are [0 



Sff(A) = SZZ(X) = 1 (27) 

sinh(A) 
sinh(A + 2ij) 



+ ~ v ; ~ +v ; sinh(A + 2i 7 ) v ; 



where the rapidities Aj(A = Ai — A2) are related to the momenta kj by 



shiKV2-j7/2) 
ex P^J" s inh(V2 + i T /2) m 
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Finally, the boundary scattering matrices are 

,(/), _ l + [P 1 ( j ) -2cos( 7 )]e^ 
1 + [PS - 2 cos(7)]e-^ 



site**®) = r:^ : : : .;. m 



Sr(k.pW) - 1 + ^ ) - 2cos ^ e ' ikj c ^^ (32) 



where 



and 



pW _ pW _ cosh(g_/2+n) 
p (l) = I M+ 1- cosh(^/2) /ggN 

' lr > p (2) = cosh(g_/2+n) p (2) = sinh(g_/2-2i 7 ) 1 j 

1+ cosh(£_/2) ' 1- sinh(5_/2-i7) 

I p(l) _ p(l) _ cosh(g + /2+i 7 ) 
p(m) _ J ~" cosh«+/2) /g 4 N 

rCT I p(2) = cosh(g + /2+i 7 ) p (2) _ sinh(g + /2-2» 7 ) 1 j 

L r + cosh(£+/2) ' r r- sinh(£ + /2-r-y) 

So far we managed to solve the charge degrees of freedom of the system but we still have to 
diagonalize the spin sector associated with the scattering matrices Sij(ki, kj), S\(ki, P®) and 
Sr(ki, Pffl)- These amplitudes, however, are easily related to those of the six-vertex model and 
the spin part of the problem is reduced to the diagonalization of an inhomogeneous 6-vertex 
model with open boundaries [18 1 . In the course of solution one has to introduce a second Bethe 



ansatz for the spin rapidities fj,j,j = 1, • • • , iV+. Since this problem has been discussed in many 

different contexts in the literature JT7L ^, [Tj| we restrict ourself to present only the final Bethe 

ansatz results. For the four possible boundary case, we find that the pseudomomenta Xj and 

the spin variables fij satisfy the following nested Bethe ansatz equations 

•2L ah 



sinh(A j /2 - ij/2) 
sinh(A i /2 + i7/2) 



„ n c v 1f sinh(Aj -fi k - ij) smh(Aj + fi k - ij) . 

f {Xj, 4±J =11 . rrr — r -r : rrr — : ;— r^r, J = t, ■ ■ • , iv e 

^ smh(Aj - Hk + ij) smh{Xj + [i k + ry) 

(35) 



-pj- sinhQj - A fc - zt) sinh(/x i + A fc - ry) _ _ ^j- sinh(/x j - /x fc - 227) sinh(/i j + - 2ry) 
^ sinh(//j — A fc + 27) sinh(/ij + X k + ij) J ' 11 sinh(//j - + 2*7) sinh(//j + - 2ry) 

k = 1 

j = l,--,N+ (36) 

where the boundary factors F(Aj,£±) and G(nj,£±) are given by 

, cosh(A J /2 + ^ 7 /2 - £-/2) co S h(A 3 -/2 + 2 7 /2 - £+/2) , , 

1 J '^ ±J cosh(Aj-/2 - i 7 /2 + ^_/2) cosh(Aj-/2 - i 7 /2 + ^+/2) 1 J 



cosh 2 (nj—i'y) 
cosh 2 (fj,j +i~f) 



sinh(^t 



sinh(|U, 
cosh (/j, 



cosh (/j, 
sinh(/ij — 



-t27+g + ) cosh(/jj— i-y) 



+«27-g + ) coshQz,- +17) 
— 17) sinh(/ij— 2«7+5-) 



I — m— 1 
I = l, m = 2 

Z = 2,771= 1 



(38) 



+17) sinh(/ij— £_+2?7) 

2 7 +g_) sinh(/jj-2i7+g+) ^ _ _ 2 



, sinh(/ij+i27— sinh(/ij+2i7— 

and, in terms of the rapidities Aj, the eigenvalues E^ l,m \L) are given (modulo additive con- 
stants) by 



N, 



E {l > m \L) = 



2 sin 2 (7) 



(39) 



fr[ cos(7) - cosh(Ai) 

We close this section commenting on two special open boundary conditions. The quantum 
group is obtained from the results for I = m = 2 in the limits £_ — > — oo and — > +00. 
We remark that there is a more transparent way to derive such boundary condition, however 
we have to use a different i?-matrix embedding. For further details see Appendix A. Other 
interesting boundary, concerning critical behaviour, is the free-boundary condition. We note 
that this case is achieved by setting = = in + in the model I — m — 1. 



4 Conclusions 

We have completed the analysis of the integrability of an interesting supersymmetric Hubbard- 
like model in the presence of boundary fields. This was accomplished by first deriving diagonal 
solutions of the reflection equation associated with a particular U q [Osp(2\2) invariant i?-matrix. 
This leads us to four boundaries conditions for the corresponding one-dimensional Hamiltonian, 
which have been diagonalized by the Bethe ansatz approach. Quantum-group invariant solu- 
tions have been discussed either as a special limit of the free-parameter £± or by the analysis 
of other possible i?-matrix embedding. 

The Bethe ansatz equations of section 3 provide us a tool to compute the thermodynamic 
behaviour and the finite-size corrections to the spectrum of the system. In principle, this 
allows us to determine the scattering of the physical excitations and the bulk and the boundary 
critical properties of the underlying field theory. These computations could be of interest as an 
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alternative way to rederive the results of ref . [ 15| for the integrable double sine-Gordon model. 
This also opens the possibility to obtain extra information concerning the operator content of 
this system which should provide further insight to the problem of tunneling in quantum wires. 

Finally, we mention that one possible generalization of this work is to investigate operator 
valued solutions of the reflection equation associated with the Osp(2\2) i?-matrix (5) f20| , |19| , 
21]. This will leads us to an electronic system with Kondo impurities |2JJ |19[ which hopefully 



could be the lattice analog of an interesting double sine-Gordon model with Kondo impurity. 
We plan to investigate these problems in future publications. 
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Appendix A : Other i?-matrix embedding 
The purpose of this appendix is to discuss an extra i?-matrix embedding for the Osp{2\2) 
vertex model. The i?-matrix has the same structure of equation (5) but with new weights 
Q-q/3(A) for a 7^ (3, namely 



q a e x + 1 



a - j3l_ 



1 + e 

a - 1- q 2 
1 + e 



1 i 1- Z 
~X + 7* — 



q - + V4<5 



1 a,/3 



where S\ = — e 4 = q, e 2 = —e 3 = i and a is defined by 



a- 



a-\ \<a<2 
a+\ 3<a<4 



a = j3 
a < (3, 

a > (3 



(A.l) 



(A.2) 
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V 



(A.3) 



This i?-matrix satisfies the properties (2,3), but now the crossing matrix V is 

/0 
iq- 1 
-iq 
\q / 

We note that some of the new Boltzmann weights a a p(X) have indeed a different functional 
form as compared to those of equation (6). For periodic boundary conditions, however, we have 
checked that such differences are not important as long as Bethe ansatz analysis is concerned. 
The corresponding Bethe ansatz equations of this "new" vertex model (or associated quantum 



spin chain) are precisely the same as that found in ref.|10|. The situation for open boundary 
conditions is, however, not so rich as in section 2. Although we managed to find two classes 
of diagonal i^-matrices solutions, none of them possess a free-parameter. The first solution is 
the standard quantum-group invariant one 



K W (X) = Id 



(A.4) 



while the second class is given by 






V o 





A 2 (\) 








A 2 {\) 




o \ 






(A.5) 



where 



A 1 (X) = (e x + eiq 3 )(e x - eiq~ 3 ) 
A 2 (X) = (e~ x + eiq 3 )(e x - eiq- 3 ) 
A 3 (X) = (e- x + eiq 3 )(e~ x - eiq~ 3 ) 



(A.6) 



where e = ±1. The Bethe ansatz solution for such open boundaries follows closely the steps of 
section 3. The only difference is concerned with the bulk scattering matrix of the spins degree 
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of freedom. Now, this matrix possesses the quantum-group invariant form: 

S++(A) = SZZ(X) = 1, (A.7) 

S+ + W = S+-(\)= g^j , (A.8) 

+ ~ v ; ~ +v ; sinh(A + 2i 7 ) v ' 

SZt(X) = e* .^- V (A - 9) 
smh(A + 2ry) 



sinh(2ry) 



Si" (A) = e~ A - /A V ,y r . (A.10) 
+ ~ v ; sinh(A + 2i 7 ) v ' 

Having this information, it is easy to derive that the Bethe ansatz equations associated with 
the first boundary (A. 4) are of quantum-group type, i.e. F(\j,£±) = G([ij,£±) = 1. Similarly, 
the Bethe ansatz equations for the second boundary (A. 5) gives us F(\j,£±) = / e _(Aj)/ e+ (Aj) 
and G(fj,j,£±) = 1 where 

! cosh[A/2-»(7+7r/4)] m 
cosh[A/2+i( 7 +7r/4)] (All) 
sinh[A/2-t( 7 +ir/4)] _-, 1 ' j 

sinh[A/2+i(7+7r/4)] 

We remark that (A. 11) can be recovered from our previous results for the "mixed" boundary 
conditions I — 1, m — 2 and I — 2, m — 1 via fine tuning of the parameters £±. The conclusions 
of this appendix suggest that such different embedding may be formulated as a twisting of 
Deguchi et al [ffT[] original solution It would be interesting to explore this possibility further 



since this may lead us to new integrable multiparametric spin chains [22 . 
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